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RADIATION, CONVECTION AND CONDUCTION COEFFICIENTS IN 
SOLAR COLLECTORS 


H. TABOR 
The National Physical Laboratory of Israel, Jerusalem 


ABSTRACT 
Heat losses from Solar collectors are usually computed on the basis of figures given by 
Hottel and Woertz!. The present paper suggests, after a study of the literature on convection, 
new Coefficients for convection heat transfer as well as drawing attention to the value of the 
total hemispherical emissivity of glass, and to the edge losses that can occur in collectors 
of finite size. 


In the design of solar energy collectors, particularly of the flat plate type, the coe- 
fficients of radiation, convection and conduction heat transfer from the heated plate 
must be known in order to compute the efficiency of collection. Most workers in 
the field have based their calculations on figures given in the classical paper of Hottel 
and Woertzl. 

The present paper suggests revised figures for these coefficients. 


BASIC EQUATION 


Hottel and Woertz have written down the basic equation of heat loss from a heated 
plate covered by one or more air-spaced sheets of glass. 


TA —T,4 
jmlmegay i) eh tal (1) 
rg ge 


Where g: = rate of heat loss per unit area of the receiver 
he = convection loss coefficient 
Ty, 7, = temperatures (absolute) of the plate and the cover glass respectively 
Ey, Eg = the emissivities (total hemispherical) of the plate and the glass cover 
respectively 
o = Stephan-Bolzmann constant. 


The first term on the right hand side of equation (1) represents the heat transferred 
by convection whilst the second term represents that transferred by radiation. A 
similar expression represents the heat transfer between two glass plates with the 
difference that Ep is replaced by Ey in the radiation term. 
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For the convection coefficient, the authors have written 
he = c (To — T,)'/4 


where c is a constant dependent upon the angle of tilt, so that equation (1) becomes: 


o(T)4 — T;4 
gi = C(T) — T,)5/4 + eile (la) 


BET 


for the transfer from the plate to the glass, and 


pep 4 
gi = ¢ (Ty — T))5I4 + ee (1b) 
pe 


for transfer from glass 1 to glass 2 and similarly for any further glass sheets. 
For the outermost sheet (temperature T,) the heat loss is 


gi = hw (Tn — Ta) + Ego (Tn4 — Tad) (Ic) 


assuming the effective radiation temperature of the sky to be approximately equal 
to Ta, the temperature of the outside air. hy is the loss coefficient due to wind which 
is known if the wind velocity is known. 

In attempting to find a suitable numerical value for the coefficient c the authors 
observed that the literature gave coefficients varying by as much as 50% as shown 
in their Figure 15. They therefore determined c by measuring the heat loss from a 
collector employing 3 cover glasses with 1’ spacings and tilted at 30° to the horizontal. 
By inserting the known values of the emissivities Eo, Eg and solving the simultaneous 
equations of the type (1a) (1b) (1c) they found a value of c = 0.166 when temperatures 
are measured in °R (Fahrenheit absolute) and the other units of measurement are 
the foot, the hour and the BTU. This value is slightly above the value suggested by 
Fishenden and Saunders? (0.147) and well below the value suggested by Wilkes3 
(0.22). 

- This experimental method of determining c involves one serious potential source 
of error: if the emissivity of the plate or the glass is not accurately known, it can 
introduce a large error in the value of c. For in equations (1a) and (1b) when a nume- 
rical value of 0.166 is substituted for c and a value of 0.172 x 10-8 BTU/hr, ft2, °R4 is 
substituted for o, it is found that the radiation term is about three times as large as 
the convection term for temperature differences of the order of 40°F*. From this 
it follows that an error of 1% in say Eg or Ep in equation (la) wil! result in an error 
in the value of c of about 3% ,whilst in the case of equation (1b) an error of 1% 
in Ey will result in an error in c of about 6%. 


* The Hottel-Woertz experimental collector had a temperature rise of 104-126° F for a system of 3 air 
gaps. 
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It is unfortunate that the authors chose a value (0.96) for the emissivity of glass E 9 
based upon the principle that the emissivity of a surface which is opaque to thermal 
radiation — as is glass — is the complement of the reflectivity, the latter being about 
4%. The last statement is probably incorrect. In the visible spectrum the normal 
reflectivity is about 4% and would appear to be about 6% in the heat spectrum, 
since the normal emissivity which is complementary to the normal reflectivity has 
been measured by Schmidt and is given as 0.94. As is well known, glass does not 
obey Lamberts’ cosine law and like most non-metals the reflectivity increases with 
angle of incidence — as predicted from the electromagnetic theory of light. As a 
consequence, the emissivity in directions other than normal to the surface is generally 
lower than the normal emissivity resulting in a total hemispherical emissivity lower 
than the normal emissivity, the value given by Schmidt4 being 0.88*. 


If this value is substituted in the formulae (la) (1b) (1c) applied to the measured 
heat loss in the Hottel and Woertz experiment, the value of c becomes 0.21 or about 
26% more than computed by the authors and, as is shown below, is probably very 
near the correct value (See Appendix II). 


Jt is interesting to note that, in many calculations of the performance of flat 
plate collectors, the use of the value 0.166 for the convection constant (for 30° tilt) 
results in practically no error if the value 0.96 is used for Ey. For in the case of a 3 
glass collector the considerable errors in the convection and radiation terms cancel 
out entirely, whilst for 2 and 1 glass systems a small error is introduced. This argument 
of cancellation of errors applies to collectors with receiver plates having emissivities 
of the order of 0.95. However, when a receiver plate of low emissivity is used, there 
is no corresponding cancellation of errors. As shown elsewhere5, selective black 
surtaces have been prepared with emissivities of the order of one tenth that of normal 
black surfaces. In such cases the radiation term of equation (la) may be only one 
third the magnitude of the convection term so that an error of 26% in the value of 
c can have a large influence on the accuracy of the total heat loss calculation. Indeed 
it was the experiments in this laboratory on a selective black collector which led to 
a study of the convection and radiation coefficients. Experiments on such a collector 
with one sheet of glass and a plate emissivity Ey) =0.11 showed an equilibrium 
temperature of 342 °F under a solar insolation of 314 BTU/hr, ft? whereas the calculated 
equilibrium temperature using the methods of Ref 6 (which employed the Hottel 
and Woertz coefficients) was 373°F**. 


CONVECTION BETWEEN PARALLEL PLANES 


From considerations of dimensional analysis, it is customary to write the natural 
convection transfer between two parallel planes as: 


* This is the value at 90°C. At room temperature the value can be expected to be very slightly higher. 


** The curves shown in Figures 5, 6 and 7 of Ref. 6 now require recomputation. 
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Nu = C (Gr.Pr)" (2) 


where Nu = Nusselt number (dimensionless) = heL/k == (qcL)/(k A T) 

he = convection heat transfer coefficient (BTU/hr, ft2, °F) (kcals/hr, m2, °C) 

L =characteristic dimension. In the case of parallel planes this is the 
distance between the planes (ft) (m) 

k = thermal conductivity of the air (BTU/hr, ft2, °F/ft) (kcals/hr, m2, °C/m) 

qe = convection transfer per unit area (BTU/hr, ft?) (kcals/hr, m2). 

T =temperature difference between planes (°F) (°C) 

Gr = Grasshoff number (dimensionless) = gy A TL3 02/u2 

Pr = Prandtl number (dimensionless) = cy pu/k 

Gr.Pr = (gy AT L3 @? cp)/(uk) 

g =acceleration due to gravity (4.17 x 108 ft/hr2) (1.27 x 108m/hr?) 

y =coefficient of volumetric expansion of air (°F-1) (°C-1) 

o = density of air (1b/ft) (kg/m3) 

Cp = Specific heat of air at constant pressure. (BTU/1b, °F) (kcals/kg, °C) 

f= ©-= viscosity of air (1b/hr.ft = 2.42 x viscosity in centipoise) (kg/hr,m = 
3.60 X viscosity in centipoise) | 

C, n = constants determined experimentally, for particular ranges of (Gr. Pr). 


Note: (1) that the product (Gr.Pr) may be written as aL3;AT where a = gyo2cp/uk— 
which is a property of the air (ft-3 °F-1) m-3 °C-1). This form is convenient for 
computation when a is tabulated*. 


: 
| 
| 
| 


Note: (2) that the Nusselt number, in this particular case of parallel planes, is the 
ratio of the convection heat transfer coefficient h- to the conduction heat transfer 
coefficient k/L across the gap if no convection were present. 

Equation (2) may be written: 


Nu = heL|K = C(Gr Pr)" = C(aL3 AT)" (3a) 


or he = C(ka”) — = (3b) 


A number of authors have reported on actual experiments carried out to determine 
the constants C and whilst others have sought to fit curves to these experimental 
results. Because of different constructions and experimental techniques the results 
do not all agree. Thus Mull and Reiher7 (1930) carried out an extensive set of expe- 
riments to which Jakob (1946) fitted two lines for the case of horizontal planes, 
heat flow upwards. 


* Some authors omit the Prandtl number from Equation (2) because it is substantially constant 
for air (0.72) over a wide range of temperature. Thus Gr = (a L3AT)/0.72) See Appendix I for tabu- 
lated values of a. 
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1. Nu =0.195 Grt/4 = 0.213 (Gr.Pr)1/4 for Gr = 104-106 (4) 
2. Nu =0.068 Grt/3 = 0.077 (Gr.Pr)1/3 for Gr > 105 (5) 
the lines crossing at about Gr = 3.7 x 105. 


For vertical planes, Jakob suggested: 
1. Nu =0.18 Gri/4 (A/L)-9 for Gr =2 x 104-2 x 105 (6) 
2. Nu =0.065 Gri/3 (H/L)-1/9 for Gr =2 x 105-107 (7) 


where H is the vertical dimension. The (H/L)-1/9 has been questioned by De Graaf 
and Van der Held!9 and Jakob!! has admitted its uncertainty. 

De Graaf and Van der Held9 carried out a number of experiments measuring the 
heat transfer at various angles. They stated that their results were very similar to 
those of Mull and Reiher but disagreed with the conclusions of Jakob. Their results 
are as follows: 


a) Horizontal planes, heat flow upwards 
1. Nu lfor.Gr < 2000 


2. Nu = 0.507 Gr°.4 for Gr = 2000 — 50,000 (10) 
3. Nu = 3.8 (almost constant) for Gr = 50,000 — 200,000 (11) 
4. Nu = 0.0426 Gr?-37 for Gr > 200,000 (12) 


b) 45° slope, heat flow upwards 
1. Nu =1 for Gr < 4000 


2. Nu = .0503 Gr-37 for Gr 4000 — 50,000 (14) 

3. Nu = .0372 Gr-37 for Gr > 50,000 (15) 
c) Vertical planes 

1. Nu =1 for Gr < 7000 (16) 

2. Nu = 0.0384 Gr0-37 for Gr 10,000 — 80,000 (17) 

3. Nu = 0.0317 Gr?-37 for Gr > 200,000 (18) 


It will be noticed that there is some discontinuity in the Nu — Gr curve in all 
three cases chosen. Wilkes and Peterson} using a gap width of 0.302 ft and plate 
temperatures of 60°F and 40°F recommended h- = 0.18 A T1/4 for vertical plates, 
0.207 A T1/4 for 45° tilt (heat flow up) and 0.256 A T1/4 for horizontal (heat flow up). 
This is equivalent to writing Nu = C (Gr.Pr)!/4 where C = 0.25, 0.287, 0.355 for 
the three cases, with the implied limitation that this form may be used only for Gr.Pr 
values around 106 the value used in the experiments. 


Fishenden and Saunders,? for vertical gaps, suggest 
Nu = 0.15 (Gr.Pr)!/4 = 0.138 Gr1/4 for (Gr.Pr) 104 — 106 (19) 


whilst for (Gr.Pr) values > 106, which generally means larger gaps, they deduce 
Nu = 0.05 (Gr.Pr)!/3 = 0.0448 Gr1/3 (20) 
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In 1954 the Housing and Home Finance Agency in Washington published Housing 
Research Paper No. 32 entitled “The Thermal Insulating Value of Airspaces’’ being 
a project that the agency sponsored at the National Bureau of Standards. This — 
paper, designed primarily for builders and architects, does not mention Grasshoff, 
Prandtl or Nusselt numbers but these were used in bringing all the experimental 
results on to a single log-log graph (their Figure 8) which shows /¢L against ATL3 
for a single mean temperature of 50°F. For convenience in computation, the paper _ 
also gives a number of graphs showing heso (i.e. at SO°F) against gap spacing L for | 
various temperature differences AT, five such graphs covering the following cases: 
horizontal planes heat flow upwards; 45° slope, flow upwards; vertical planes, 
flow horizontal; 45° slope, flow downwards; horizontal planes, flow downwards; 
These five cases are also shown on their Figure 8. It is possible to fit straight lines 
to the curves of Figure 8, the fits being very close for‘the horizontal upwards and 45° 
upwards cases, much poorer for the vertical and 45° downward cases and not possible | 
for the horizontal downwards case. 


Thus we find: 
d) Horizontal planes, heat flow upwards | 
heso = eee A T°.281 for Gr f 104 — 107 21 | 
050 = 79.187 or Gr from rs (21) 


where / is in inches but he is in BTU/hr ft2. 


e) 45° planes, heat flow upwards. 


heso = a A T9-310 for Gr 104 - 107 (22) | 

f) Vertical planes. | 
(i) Aesop = ae A T°9.327 for Gr 1.5 X 105- 107 (23) | 

(it) Aesq = 0.0917 A 79.381 19.143 for Gr 1.5 x 104-1.5 x 105 (24) | 

| 


By introducing the conductivity of air at 50°F and arranging the terms in aon- 
dimensional form, these results may be written in generalised form without being 
limited either to a single temperature, or to English units. Thus (d) becomes: | 


Nu = 0.168 (Gr.Pr) 0.281 = 0,152 Gr 0.281 (25) | 


e) becomes: | 
Nu = 0.102 (Gr.Pr) 0310 = 0,0925 Gr 0.310 (26) | 


fi) becomes: 
Nu = 0.0685 (Gr.Pr) 9.327 = 0.0616 Gr 0.327 (27) 


and (fii) becomes: 
Nu = 0.0369 (Gr.Pr) 0.381 — 0.0326 Gr 0.381 (28) 
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An examination of all these figures shows a moderate correspondence between 
the results of Mull and Reiher as expressed by Jakob, of De Graaf and Van der 
Held and of Report 32, with the latter taking an intermediate position for the case 
of horizontal planes. Table I shows values of Nusselt numbers computed from the 
various formulae, for standard Grasshoff numbers Gr. The value Gr = 3.86 x 104 
has been inserted as this corresponds approximately to the conditions existing in 
the Hottel and Woetz (H and W) experiment described earlier (mean temperature 
Tm ~ 100°F, L = 1/12 ft, AT ~ 40°F per gap). The H and W values of the Nusselt 
number for various angles of tilt are computed from their Figure 15 and have been 
limited to the range Gr = 104 — 105. It is seen that the Wilkes and Peterson results, 
which are for a single value of Gr approx. 1.39 x 106 are much higher that the other 
results, which Jakob attributes to the error of the conducting studs and the absence 
of a guard ring. As Jakob points out, the work of Mull and Reiher gives results 
io be expected in the case of an enclosed partition cut out in the middle of an extented 
air layer between two parallel surfaces, with a practically non-conducting border. 
The same comment can be applied to the later work of Report 32 where a guard 
ring is used. Clearly experiments carried out on finite size plates without a guard 
ring — such as Wilkes and Peterson — can be expected to give higher conductance 
values. 

The present writer feels that, without detracting from the work of other experi- 
menters, the results given in Report 32 are probably the most reliable and consistant 
set when applied to fairly large plates. A correctlon, at present unknown, needs to 
be applied when dealing with small plates, The writer’s preference for this work is 
based upon the fact that some 146 tests were carried out with gaps varying from 
about 4” to 3.4” (A/L 18-96) and with both low and high emissivity surfaces. The 
spread of the results on Figure 8 of that report and shown in the present Figure 1. 
is surprisingly small except for the case of horizontal gaps, heat flow downwards, 
where unstable conditions are to be expected. 

Another satisfactory feature is that, for large Gr numbers (which generally means 
large gap spacings L) the index n in equation (2) is found to approach the value § but 
does not exceed it. This appears logical since from equation (3) the heat flow coe- 
fficient becomes independent of L for n = 4 as is to be expected if one considers the 
heat transfer in two stages — from the hot plane to the air and from the air to the 
cold plane. Jakob, Fishenden and Saunders and others have recognised this, though 
De Graaf and Van der Held give values of 7 as high as 0.37. 

A further point to be remembered is that for convection experiments, edge effects 
can be quite important and it is thus unlikely that the results obtained in a real solar 
collector will be identical to those obtained in a careful laboratory set-up. The aim 
of convection computations is, therefore, to obtain as consistant resu'ts as possible 
even if the absolute values are somewhat different, for it is then possible to estimate 
quite accurately the effect of a change in design parameters, thereby facilitating the 
approach to an optimum design. 
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TABLE I 


Nusselt numbers A-L/k for parallel air gaps 


Gr. ——> 103 = 104 3.86 x 104 105 106 107 
6.15* 
Horizontal Jakob(MandR) 1.0 1.95 23 3.46 68 14.6 
Heat flow Gand H 1.0 2.02 3.46 3.8 7.1 16.6 
upwards Report 32 1,0. t52:02 2.94 3.85 E35 14.0 
W and P — — — — 10.3 
Hand W — 1.81 2.54 Sle — 
45° slope, 
heat flow G and H Te Oo mealeo 2.49 2.64 6.17 
upwards Report 32 Ov 16! 2.44 3.28 6.78 
W and P —_ — — — 8.35 
H and W — 1.46 2.06 2.60 
Vertical 
heat flow Jakob (M and R)** 1.81 2.29 4.65 
horizontal Gand H 1.0 1.16 1.84 — 5.26 
Report 32 COM S085 * ss 1283 2.62 eo, 14.1 
W and P Wes 
FandS HRO) ees 1.91 2.45 4.35* 9.65 
H and W 1.26 1.79 2.24 4.48 
MandR Mull and Reiher * Transition values 
GandH De Gratf and Van der Held ** Assuming A/L = 20 
WandP Wilkes and Peterson *** Poor fit of equation f (ii) at this point. 


Hand W Hottel and Woertz 
FandS _ Fishenden and Saunders 


Figure 8 of Report No. 32 has been redrawn in Figure 1 with the addition 
of Nu and Gr numbers. Furthermore, by introducing numerical factors F, and F2 
which allow for the variation with mean temperature of Gr and Nu numbers (for air), 
the graphs can be used over a wide range of temperatures. 

The graphs, when the ordinates and abscissae are given as non-dimensional Nu 
and Gr numbers are perfectly general. However, for convenience in use the ordinate 
is also given as h,/* iastead of h-L/k and the abscissae as A773 instead of aATL3 
(See equation 3a). These are now dimensional and their use is permissable for single 
values of k and a, i.e. at one mean temperature, which is 50°F in Report 32. 

For hA¢l/k, at any temperature T, we can write Fyh¢l/ks59 
where Fy = kso/k, and for a-ATI3 we can write F\a5,A4713 
where F; = a;/asy We can thus use the dimensional scales of hl and AT/3 for any 
chosen temperature provided that we insert the values of F, and F, for that tempe- 
rature. 


The method for using Figure 1 is given in the example below the figure. 


* The small / is a departure from consistant units (i.e. inches instead of feet and cms instead of 
meters) because of the usual gap spacings being of the order of one inch. 
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EFFECT OF MEAN TEMPERATURE 


The influence of mean temperature on the values of the convection coefficients is 
rather small in most practical cases. Equation (3b) shows that the coefficient is pro- 
portional to (ka”) when the Nusselt number is a single power function of the Grasshoff 
number, i.e. on the linear portions of the curves of Figure 1. As indicated in equations 
(25), (26), (27), m has the values 0.281, 0.310, 0.327 for the cases of horizontal planes 
(heat flow up:) 45° tilt planes (heat flow up:) and vertical planes respectively. A 
check on the variation of the quantity ka” has been made in Appendix I, and shows 
that for these three cases, the coefficient decreases by 0.6-0.8 % for each 10°F increase 
in mean temperature. For heat flow downwards the effect must be checked by use 
of Figure 1 at several mean temperatures at the required A7/3 values: It will be 
found that for 45° slope (downwards) the effect of mean temperature is rather small, 
whilst for the horizontal case (downwards) the coefficient increases with rising mean 
temperature by roughtly the same order of magnitude as the variation for heat 
flow upwards. 


EFFECT OF GAP DISTANCE 


The effect of changing the gap distance / can be obtained, for any set of conditions, 
from Figure | by carrying out the computation, as given in the example, for several / 
values. In the case of values falling on the linear portions of the graphs the effect 
can be determined at once from equation (3b), i.e. Ae is proportional to L3”-1. Thus, 
for the cases of horizontal planes (heat flow upwards:) 45° tilt, (heat flow upwards:) 
and vertical planes, where n has the values 0.281, 0.310, 0.327 respectively, the 
influence of L is L-9-157, 1-09.07, [-0.019 or increasing L by a factor of two reduces the 
convection coefficient by 10%, 5% and 1.3% respectively. The use of large gaps, with 
their attendant increased shading losses and increased box costs is not justified 
except for very efficient collectors. Thus for collectors in the tilt range 30-60°, a 1” gap 
as used by Hottel and Woertz is very reasonable. For collectors using selective 
black receiver plates, where the major heat loss is by convection, a slightly greater gap 
spacing may be justified. (For heat flow downwards, at 45° tilt, increasing the gap 
distance has sometimes a positive and sometimes a negative effect on the convection 
coefficient whilst for the perfectly horizontal case, heat flow downwards, this coe- 
fficient is greatly reduced because most of the heat is transferred by pure conduction). 
For the designer working continuously with one gap dimension, and near a given 
temperature, it is possible to prepare specific graphs from Figure 1. 


Thus Figure 2 has been prepared from Figure 1 for the case of 1’ (2.5 cm) air gap 
and a mean temperature of 100°F (38°C) with curves for 30° and 60° slope added 
by interpolation. For mean temperatures other than 100°F or for a 2” gap instead 
of 1” gap, the approximate correction factors can be used as shown. The heat transfer 
downwards at 45° slope and horizontal has been added to the graph for use in cases 
where reflective insulation is considered for the rear of flat plate solar collectors. 
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Generalised graph of heat flow by convection between parallel planes. (Derived from Figure 8 of Ref. 14). 


Ex. Determine convection coefficient for parallel planes tilted 45° (heat flow upwards) gap spacing 2 ” mean temperature 150°F temperature 


difference 100. From Table, F; = 0.45, F, = 0.850. 
“FF, ATI = 360. Foe = 0.98 she = 0.577 BTU/hr, ft2,°F 
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Generalised graph of heat flow by convection between parallel planes. (Derived from Figure 8 of Ref. 14). 
Ex. Determine convection coefficient for parallel planes tilted 45° (heat flow upwards) gap spacing 2’’, mean temperature 150°F temperature 
difference 100. From Table, F; = 0.45, F, = 0.850. 
.. FATB =360. Fohce=0.98 .*.he = 0.577 BTU/hr, ft?, °F 
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Where the designer is working with one slope and one gap size it is more convenient 
to prepare curves showing the actual heat transfer ge =he AT against AT rather 
than h-. On the same graph it is possible to plot the radiation heat transfer if the 
effective emissivity between the planes is known. If qc is plotted as the positive ordinate 
and qr is plotted as the negative ordinate, the sum of the ordinates gives the total 
heat transfer between the planes. Figure 3 is an example of such a set of curves, in 
which the radiation heat transfer has been plotted for two cases — one, the case of a 
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receiver plate, of emissivity 0.1 with a “black” (or nearly black) cover, (ii) the case 
of two sheets of glass assuming Ey = 0.88 and FE. = 0.786. 

As indicated in Ref. 6 this type of graph is particularly useful for computing 
performance of solar collectors having multiple glass panes. 


III. REAR CONDUCTION LOSS 
The heat loss from the rear side of a solar receiver plate is calculated by the usual 
methods 


Ix = UR(Tp — Ta) (29) 
Where Up = overall rear loss heat flow coefficient 
1 1 (30) 
1 1 d 1 
Upset a = teat 
: hp Agr k Agr 


where d is the thickness of insulation of conductivity k, hoz is the outside loss coe- 
fficient, and Ap is the rear conductance through the insulation alone. The value of 
Nor is usually not known very accurately but will generally be about 2-4 BTU/hr, 
ft2, °F (5-10 kcals/hr, m2, °C) and may be even less for a solar roof construction 
where there is no wind at all on the rear side of the collector, However, as /p is 
usually much smaller than /gp a moderate error in the value of the latter coefficient 
is not very significant. 

The purpose of the present note is to draw attention to a point generally ignored 
in the literature and which can result in a serious error, or a waste of insulating 
material, in the case of collectors of small to moderate dimensions. 

The formula (30) is that applying to a parallel slab of insulation of a finite thickness 
but infinite extent. When the slab is not infinite, the heat flow lines are curved at 
the edges resulting in a greater heat flow per unit area than for an infinite plate. 

Thus, the special case of a finite thin rectangular plate of dimensions L x M and 
surrounded on both sides and at the edges by a rectangular layer of insulation of 
thickness 6 has been treated by Langmuir, Adams and Meikle (13) and is quoted 
by McAdams (12) and in Ref. 6. 

Thus the heat flow from the plate at temperature Ty to the outer face of the insu- 
lation assumed to be at uniform Fait T, is . . 


in)! 
(L) Gi ) oe 
or the effective area for computing heat flow from the two sides of the plate is inc- 
reased by the “edge correction factor” F, where 


Q =2LM [1+0.93(+ a a) + 0.175 ae (31) 


8b waigd 56 6 
Fs SU 000. 4) Ol ae (32) 


The second term represents the effect of the edges, the third term the effect of the 
corners. The latter term is, for all usual practical cases quite negligeable. 
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Equation (32) may be written, omitting the effect of the corners as 


Pas ps 


De 
where L is the “characteristic” size of the plate and is given by 
1 1 1 
~=F5taZ (33) 
— L 
Z M 
and S has the value 0.93 for this case. The term S6/L may be called the “edge 
addition factor” and the constant S$ the “‘edge addition coefficient”. 


Thus for a plate 3 ft x 3 ft with an insulating layer all round of say 4” thickness, | 


the edge correction factor F, is 1.207 or the heat flow per unit area is 20.7% * more | 


than for an infinite plate. The question arises, what happens if the edge insulation 
is not as thick as the rear insulation, a condition found in many practical designs 
where the rear insulation is of the order of 2-4” of insulating material whereas the 
edge insulation may be a wooden frame of perhaps only 1” thickness. Furthermore, 
in the special case chosen by Langmuir the conductance on both sides of the plate 
is equal, a condition not necessarily applying in a collector system. 

A rigorous mathematical analysis of the generalised case of a real collector is 
extremely difficult though it is possible for a hypothetical collector: i.e., for a collector 


in which the heat flow from the top face of the heated plate is also purely by con- 


duction. 


The total heat flow from a plate with insulation thickness 6’ on one side (the “top’’) | 


and 6 on the other side (the “rear’’) may be written . 
1 1 ) 
Q=IM(T, —-Tk(G + FL +S] (34) 
7 


where the square bracket is the edge correction factor for a plate of finite size and 
the edge addition coefficient S may be a function both of the ratio 6’/d and the ratio 
b/6 where b is the thickness of the insulation at the edges. See Figure 4. 


PLATE 


INSULATION 


(a) LANGMUIR CASE (6) GENERAL CASE 
Figure 4 


* If the heat flow on one side had been taken as LM k/AT the flow on the other side would have 
to be taken as 1.414 k/AT, i.e. increased by 41.4% to give the correct total heat flow. 


——— 


| 
(i 
| 
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Using the theory of conjugate functions applied to the Fourrier-Laplace heat 
flow equations, Mr. H. Weinberger of this laboratory has computed the edge addition 
coefficient S for varying values of b/d for a few values of 6’/6. Note that if the heat 
flow from the “‘top”’ (thickness 6’) had been computed as if the plate were infinite in size, 
as is commonly done, i.e. Orop = LMATk (1/6’) then the loss from the other side 
and edges, which we call the total rear and edge loss is 


TABLE I 
Edge addition coefficients S, S’ 


0/6 =2 b/6 = 0 vs 1} 1 4 1/4 1/8 
S = 1.22 |e) a os Ras ais Bp ea Hy a of 0 
S’ = 1.83 1.85) TORpOD PRS Oo Sp 'op SOL AiG) 5.3 

6/5 =1 S = 0.882 0.882 0.893 0.938* 1.179 1.65 2.35 
S’ = 1.76 Lit -1,fo. . -1a5- 2.36. 333 4.7 

6/5 = 4 S= 0.611 0.615 0.616 0.62 0.71 096 1.33 
S’ = 1.83 L £4 see MOS e186) 4 813 RE a AD 


© rene = Oe 8 hs 
=LMATK[(— + +) (1+=)--] 
) om L o 


ty) 


=rmaT <i+s (1+ 2) *] (35) 
6 6 i 
= LMAT* [1 whl] (36) 
g i 


where the square bracket gives the total edge correction factor for both the top and 
bottom faces of the plate. The value of S! has also been included in Table II. 


Several conclusions can be drawn from Table II. Firstly, the edge addition coe- 
ficients S and S! are hardly affected by making the insulation thickness b at the 
edges greater than the insulation thickness 6 at the rear. But as soon as the edge 


| thickness falls below the rear thickness, the edge addition coefficients rise rapidly. 
| This suggests that when computing an optimum design, taking into account the 
- cost of increasing the edge insulation, the optimum will probably come out at near 


* Note the excellent agreement of this value (when b = 6 = 0’) with that given by Lang- 
muir (0.93) 
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the value b = 6 and this may be taken as a satisfactory starting point in any design. 
Secondly, in the region b/6 ~ 1 the influence on the edge addition coefficients of 
the ratio of the conductances on the two sides of the plate (the 6’/6 ratio) is not very 
great so that an accurate estimate of the conductance on the top side of the plate 
is not essential. This suggests that the results shown in Table IJ computed for a 
hypothetical collector, should apply, at least in order of magnitude, to a real collector. 


Thirdly, for collectors of a few feet in size employing a few inches of rear insulation, 
the edge addition coefficients are not negligeable, particularly if the b/6 ratio is less 
than 1 as often occurs. 

Thus for a collector of 6 ft x 3 ft (ie. L = 2 ft) employing 4” of rear insulation, 
the coefficient S! (taking 6/6 ~ 1) is 1.88 for b = 6, ie. for an edge insulation 
also 4’. Thus the edge correction factor (1 + S’ 6/L) is 1.31 and the 4” of insulation 
behaves as though it were only 3”. If the edge insulation was 1’’ of wood, equivalent 
to about 4” of good insulation material, the coefficient S! is 4.7 and the edge correc- 
tion factor is 1.785: the 4” of rear insulation behaves as though it were only 22”. 
(If with this edge insulation only 3” of rear insulation had been used the correction : 
factor is approx. 1.5 i.e. the rear insulation behaves as though it were only 2”: i.e. in 
increasing the rear insulation from 3” to 4” the effective increase is only from 2” 
to 24’. This stresses the suggestion already made namely to design the edge insulation 
about as great as the rear insulation.) 


WALL RADIATION LOSS 


In real collectors, energy is radiated towards the colder glass by the insulating wall 
that separates the collector plate from the glass cover. In the usual design, this, 
energy is provided almost entirely by the radiation from the hotter collector plate: 
to the wall, since the wall has a mean temperature intermediate between that of the; 
plate and that of the glass .However, if the collector plate is of low emissivity, there; 
is very little radiation transfer from the plate to the wall, so that if the wall radiates} 
to the glass cover it must “‘stea!”’ this energy from its surroundings, principally by 
conducting heat away from the plate at the line of contact with the plate. This term, 
negligeable in large collectors, becomes important for collectors less than 1 foot 
wide employing low emissivity collector plates. 


| 

No method is offered here for calculating this term. Whilst the radiation transfer| 
can be calculated for given geometrical conditions it is extremely difficult to postulate 
the temperature distribution on the wall, which is certainly not linear (between the! 
temperatures of the plate and the glass cover) due to the interaction of radiation, 
conduction and convection effects. If this distribution is known e.g. by measurement, | 
and the radiation loss to the glass is determined, it still does not tell us what part! 
has been provided by conduction away from the plate. To obtain a picture of orders| 
of magnitude, a small experimental electrically heated collector was built, the plate; 
(unpainted aluminium) being 25 x 90 cm and insulated with 5 cm of rockwood at 
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the rear and edges with 5 cm gap to the glass cover. The wall between the plate and 
the cover glass was tapered out at 45° (in anticipation of the collector being used to 
receive a convergent beam of concentrated sunlight). The heat losses were determined 
in two horizontal positions, glass facing upwards and glass facing downwards. 


The results were as follows: 


in Upwards Downwards 
Ist Plate Temp. T,)°F 340 340 
| Glass Temp. Ty °F 150 140 
| Box Temp. To°F 120 105 
in| Measured Heat Loss BTU/hr, ft? 232 143 
sal 
te Upwards Downwards 
“} Computed Radiation Loss (E = 0.05) 23 24 
| Computed Convection Loss 140 29 
uf Computed Rear Loss 58 62 
ti —— — 
_) 221 115 
idl Unaccounted Loss 11 28 


Two thermocouples embedded in the wall allowed an approximate picture of the 
) temperature profile to be obtained, from which it was computed that the wall radiated 
| about 40-50 BTU/hr (per sq. ft of collector plate) towards the glass. Some of this 
alt heat is provided by the air, by the radiation from the plate and by the edge correction 
i} in the insulation, whilst the remainder is provided by additional loss from the edges 
i of the collector plate. Thus the unaccounted losses in the computation seem to 
|) correspond — in order of magnitude at least — with the extra losses caused by wall 
4“? radiation. (That the loss is less in the upward position than in the downward may 
wf be due to: 
| The greater convection contributing a larger share of the energy, 

The actual convection (particularly upwards) may be less than the computed 

value because of the small dimensions, 

Experimental error: the plate and the box are not isothermal surfaces as presup- 

posed in the computations.) 
If the collector had been twice as wide, the effect of the wall radiation would be 
i} approximately halved: if the emissivity of the plate were greater than 0.05 as in this 
ih experiment, the wall radiation loss would again be reduced. If the connecting wall 
| had not been tapered out at 45° but had been perpendicular to the plate and glass, 
|) the loss term would again have been less due to the wall “seeing” more of the plate. 
‘|| Tt is thus seen that for most designs the wall radiation loss will be negligeable, 
'' but in the special case of narrow low emissivity collectors the term may need to 


be taken into account. 
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SUMMARISING 


a) For radiation transfer use, for effective emissivity Ee between two planes 


1 7 E, + E> 
r Jae-bal Eb + £&; — £, 


taking care to use total hemispherical values for E, and E, whenever possible. The 
value of 0.88 is recommended for glass, so that the effective emissivity between two. 
glass sheets is 0.786. The radiation heat transfer between parallel planes is 


where 77» is the mean temperature 


Taxed 
Gr = 0E¢ (T\4 — Tp4) = 40E¢ Tm? AT [1 seieseel ] 
m 


Dieh 25 


> and AT = T; — T> 


In English units, temperature in °R, 


gr (BTU/hr. ft2) = 0.00688 Ee (— 5a ™ AT [I + — -) ] 


In metric units, temperature in °C, 


gr (kcals/hr.m2) = 0.195 Ee ( Kany AT: [Licki( nat yy ] 


The square bracket can usually be ignored. 


b) 
c) 


d) 


For convection transfer use Figure 1 or graphs derived therefrom. 


For the rear loss, the edge insulation thickness b reduces the effective thickness 6 
of the rear insulation by a factor (1 + $’6/L) where L is the characteristic size 
of the plate (ZL = (LM)((L + M)) and S’ is given in Table II. For the recom- 
mended case where 6 = 6, the value of S’ is taken as 1.88. 

Thus the rear loss is: 


k aid 
L 


Vrear 


For small collectors—of the order of 1 ft wide or less — and employing low 
emissivity collector plates, an additional loss term arises from radiation from the 
boundary wall between the collector plate and the cover. 


a 


: 
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APPENDIX I 


Variation of Gr, Nu and ka" with mean temperatures 
a ert ll eee, 2 a ns eee eta 
Mean AG 10 38 65 93 121 149 Note 
Temp.T °F 50 100 150 200 250 300 


k 0.0145 0.0158 0.0171 0.0182 0.0192 0.0204 BTU hr-! ft-1°F-1 
10-14 1s 12" 0.83 O58 ode "03 ft-3 °F-1 
Fi=arlas 10 0.67 045 032 0.23 017 


| F,=kg/kr 1.0 0.917 0.850 0.797 0.755 0.711 


kruy0.281 0.834 0.811 0.785 0.760 0.730 0.710 Horizontal planes 
' Ratio* 1.00 0.974 0.942 0.912 0.875 0.853 heat flow up. Gr. = 104— 107 
krar9.310 126298122125 SIM68" 1511/5 1.06 1.02 45° Tilt, heat flow 
| Ratio** 1.00 0.960 0.925 0.886 0.840 0.808 up. Gr. = 104— 107 
krary0.327 1EGOS S36 A723 93eee 1328 1.27 Vertical, heat flow 
Ratio*** 1.00 0.955 0.916 0.867 0.827 0.790 horizontal. Gr. = 1.5 x 105 
— 107 


* (kr ar 0.281)/(k59a590-281) 
** (Kk 7a 79-310)/(k594590-310) 
*** (k 7a 7 0.327)/(K595y0-327) 


APPENDIX II 


COLLECTOR EXPERIMENTS 


A. A flat plate collector was built using a selective black receiver plate, emissivity 
E = 0.11, absorbtivity to sunlight a = 0.92. One cover glass of transmissivity B = 
0.84 and approx. the equivalent of 3 inches rear insulation of rockwool. Tilt 30°. 

The collector was exposed to bright sunshine with a peak intensity of 314 BTU/hr, 
ft2 and the equilibrium temperature T of the receiver plate was found to be 342°F. 
Outdoor temperature Tg was 80°F and there was a light wind. 


CALCULATED TEMPERATURE CF PLATE 
Of the solar radiation entering the collector, the amount reaching the plate and 
absorbed thereby is: (Ref. 6) 


qo = a8 OQ; = 0.92. 0.84. 314 = 242 BTU/hr, ft2 


which is equal to the heat loss gL(T) at equilibrium (at temperature 7). Approximately 
25 BTU are lost from the rear, so energy lost by radiation and convection from the 
plate to the cover glass is approx. 217 BTU. Heat lost from glass = 217 plus solar 
radiation absorbed by the glass, ic. 217 + 0.08 x 314 = 242. Taking the outside 
coefficient fy as approximately 4-5 BTU/hr,ft?, the glass temperature can be expected 
to be about 55°F above ambient, i.e. Tg = 135°F. 

Using the original formula (21) of Ref. 6 


T4 TA te. 
Vr) —— Qr+ Qe+ rear = 0.16AT 5/4+0.171 Eel 60° am 700. — qr —s Ta) 


| 
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and substituting the appropriate values for Ee, k, d, we find AT = 238°F 


giving de = 149 BTU/hr, ft2 
dr c= 67 
rear = 26 
Total qu = 242 


ie. Plate temperature = 373°F 


Recalculating, using the convection coefficients proposed in this report, we fiad | 
Al = 212°r 


giving Ge = 159 
Ase Sag | 
Grae = 24 | 
Total ge — 242 | 


i.e. Plate temperature = 348°F 


This is in very good agreement with the measured value of 342°F. Better agreement 
would be fortuitous siace the difference could occur by, for 2xample, a 3% error in 
measuring the solar intensity, or by taking a slighly different value for the outside 
coefficient No. 


B. The Hottel and Woertz experiments. In Ref. 1 the authors give two results of | 
their experiments on a three gap collector from which the value of the convection 
coefficient c is deduced in the equations of heat loss upwards: 


Qury = (T, — T2)5!4 + aEe (T\4 — TA) (HW11) 


or in their compounded equation (13) | 


> T — Ta a (T4 — Tq4) | 
qr) = Fae uy 1 + 1 4 n Ms, i sues (HW13) | 
IY el Se E se pappemneitce 
n+f 


where n = number of sheets of glass 
f =ratio of thermal resistance of the outer plate to that of an average inner 
plate = 0.36 for value of hg = 4.07 
E =plate emissivity = 0.95 
Eg = glass emissivity = 0.96 
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The results were: 


T°F T° F BTU/hr, ft2 % 
Exp. a 163 Bilao 57.9 0.164 
0.166 mean 
Exp. b 145.2 41.4 46.2 0.169 


Values of c were computed from the simu'taneous equations of the type (HW11) 
and substantially similar results would have been obtained using equation (HW13). 
Thus if we take equation HW13 and substitute the value c = 0.166 we get for Exp. 
(a) a heat loss g; = 57.4 and for Exp. (b) a heat loss of 46.5 BTU/hr, ft2. 


We now consider how the results would emerge if the new convection coefficients 
were used. 


In Exp. (a), the mean temperature for the three gaps is 100°F and the average 
temperature difference AT across each gap is 37.2°F.* According to Figure 2 the 
convection transfer is 0.482AT which could be written as 0.196 AT 1/4 for AT = 37.2°F. 


In Exp. (b), the mean temperature is 98 °F and the mean gap difference AT = 30.8°F 
which, according to Figure 2 gives a convection transfer of 0.4554T which may 
be written as 0.194 AT!/4, 


If now values of c=0.196 and 0.194 are substituted in equation (HW13) we obtain 
the previous loss values of 57.4 and 46.5 BTU/br,ft? if we choos: for Eg the value 
0.891 in Exp. (a) and 0.919 in Exp. (b).** 


In other words if the emissivity of glass at a mean temperature of 100°F were taken 
as 0.905 (instead of 0.96 as formerly), the convection heat losses derived from these 
experimenis would correspond almost exactly to the values suggested in the present 
report. As Schmidt’s value fo1 total hemispherical emissivity for glass is 0.88 at 90°C 
(194°F) and a slightly higher value can be expected at lower t2mperatures, these 
experiments would appear to constitute very good confirmation of the proposed 
convection values (at least in this range). 
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UTILIZATION OF DOLOMITE AS FUEL OIL ADDITIVE 


A. C. FRIEDLAND 
The Palestine Electric Corporation, Ltd., Haifa 


ABSTRACT 


The deposits of slag and ash in large oil fired boilers have a low melting point and are 
corrosive to iron and steel. Dolomite, of which a suitable grade was available locally in 
unlimited quantities, was injected into the fire boxes of some units. 

The effects of the treatment and the results of analysis of boiler deposit are described. 
Part of the dolomite reacts with free sulphur trioxide, while the unreacted raaterial decom- 
poses into calcium and magnesium oxides, which ensure protection from corrosion, The 
| resulting ash is of higher melting point and does not adhere strongly to the furnace walls 
and tubes. 


INTRODUCTION 


The accumulation of slag and ash deposits on boiler heating surfaces causes troubles 
| which may be of both a physical and a chemical nature: 


1. Fluxing (sintering) of the furnace walls. 


2. Formation of hard deposits which can not be removed satistactorily by air or 
' water lancing. 


3. Corrosion of boiler tubes, air preheaters, hangers, etc. 

The cause of the deposition of slag lies in the quantity and in the composition 
of the oil ash. Ash contents as high as 0.37 % were found and 0.1-0.2 % ash is common: 
The main constituents of the ash are: 


Sodium sulphate — Melting point 880°C 
Vanadium pentoxide — Melting point 675°C 
Ferric oxide 


with other components such as alumina, lime, silica, nickel oxide, etc. Some typical 
residual fuel oil ash analyses are given in Table I. 


Sodium sulphate and vanadium pentoxide are highly corrosive at elevated tempe- 
ratures and many tests have shown that none of the commercially available alloys 
used for the production of boiler tubes is immune from atatck when these mixtures 
are molten, that is above 670°C 1.2.3. Several research laboratories have during the 
last twenty years, studied the best ways of preventing slag deposition and ash corrosion. 
An improvement can be brought about by additives which either affect the physical 
nature of the ash (reducing its adhesive power and increasing its low melting poit) 
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TABLE I 
Composition of Fuel Lil Ash 
(1952) 
a ee ee ees ee eee ee ee 
Haifa P. ST. Reading P. ST. Reported by Babcock and Willcox (1) 
SO; 3333 44.7 37 — 67 
Na,O 4.7 15.3 10 — 25 
K,0O 15 1.8 0.1— 2 
V,0;5 36.7 4.8 1 — 35 
CaO 6.5 5.8 1 —10 
MgO 1.8 oa) 2 —3 
Al,O; 0.7 1.8 — 
Fe,0; 3.9 13.1 6 —8 
SiO, 7.0 Shs) 5 —9 
NiO 2.4 ii — 


or chemically affect the vanadium pentoxide and the sodium sulphate to produce | 
relatively inert substances. 

Several materials had been tried in the United States and Great Britain, e. g. soda 
ash, silica, vermiculite, aluminium oxide, calcium carbonate, lime, magnesium | 
oxide, dolomite, zinc oxide 14,5, 


EXPERIMENTAL | 
Since May 1955, dolomite has been added to the air stream of a number of high 
pressure boilers. This material obtained from the Havasa mine near Haifa, was 
ground by a local manufacturer* according to a prescription of the Babcock and | 
Wilcox Company so that at least 997% passed a 325 mesh ASTM sieve (average — 
particle size 0.44 microns). The recommended dolomite: ash ratio was 1:1. 
The analysis of the local dolomite is shown in Table II. 


TABLE I 
Composition of local dolomite 


Apparent Specific Gravity 1.100 
Moisture 0.1— 0.3% 
Calcium carbonate (CO;Ca) 60 — 64% 
Magnesium carbonate (CO;Mg) 27 — 36% 
Ferric oxide (Fe,03) 0.5% 

Silica (SiO,) 0.1—0.3% 
Alumina (A1,03) 0.1—0.3% 


A similar unit was run for a year without adding dolomite. The boilers were closed 
several times for inspection and on these occasions samples of fuel oil deposits were 
taken from several points and analysed. Some of the results are described in Table III. 

The first fact observed was that more ash deposited in the dolomite treated boilers 
but in a powdery, unsintered form and being soft and granular were easy to remove. 
Without additive, the deposits were of a fused, glass-like, hard and adhesive nature. 


* Messrs. Polychrom Ltd., Haifa. 
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From the analysis of the various deposits, the following facts appear to be 
significant : 

1. The total sulphur in the ash tended to be lower when dolomite is added. Much 
of the sulphur reacts with the dolomite producing calcium and magnesium sulphates 
which are infusible and non corrosive at boiler temperatures. 

2. The content of alkaline oxides was lower when dolomite was added. 

3. Vanadium pentoxide concentrated mainly in the firebox deposits. 

4. The amount of ferric oxide was very much lower when dolomite was added, 
showing that corrosion of boiler tubes and other steel parts decreased when this 
material was present. 

CONCLUSIONS 


The addition of dolomite rock would appear to be justified by several useful effects 
of which the most important are: 

1. The reaction with the free sulphur trioxide (SO3) in the flue gases to form 
calcium and magnesium sulphates which are inert and infusible. When dolomite 
is not present, SO3 reacts forming iron and sodium sulphates causing corrosion of the 
boiler parts. 


2. Unreacted dolomite decomposes into magnesium oxide (MgO) and calcium — 


oxide (CaO) at boiler temperatures, liberating CO, which escapes resulting in porosity 
of the deposits. These oxides are easily removable by water and being alkaline they 
ensure protection of the boiler surfaces from corrosion. 


3. Dolomite dilutes the ash lowering the Vanaditum Pentoxide content although 


it is not known if compounds of vanadium with calcium or magnesium are formed. 
The resulting ash has a melting point above maximum temperature of the boiler 


surfaces. 


4. Some reaction normally takes place between the boiler deposits and the re- 


fractory lining of the walls and the bottoms of the boilers, where fused material 
has been found. 

Dolomite seems to modify favourably the hardness and the sulphate content 
of the fuse product and a further study of this effect is being undertaken. 
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PAPYRUS AS A RAW MATERIAL FOR THE PRODUCTION OF 
INSULATION AND HARDBOARD* 


MENACHEM LEWIN AND AVRAHAM LENGYEL 
Institute for Fibres and Forest Products Research, Jerusalem 


ABSTRACT 


The possible utilisation of Israeli papyrus growing in the Huleh region of the Upper Galilee 
for the production of insulation boards and hardboards was investigated. Some of the 
cooking variables, the degree of refining of the pulp and the pressure employed in forming 
the boards were investigated with respect to their influence on the properties of the insulation 
boards obtained. Conditions were found in which insulating boards with properties 
conforming to British Standard specifications could be produced with a yield of about 80%. 
The laboratory results were confirmed by a full-scale mill trial, carried out at Sefen Ltd. 
in the Jordan Valley, Israel. The relatively high flexibility of the papyrus boards compared 
with boards prepared from eucalyptus was found to be related to differences in the fibre- 
length distribution of the respective raw materials and to the screen analysis of the boards. 

The preparation of hardboard by hot pressing of wet board mats was investigated and 
a temperature, pressure and pressing cycle was found which yielded boards of the standard 
hardboard quality. The influence of sizing and tempering treatments on the strength and 
water absorption of the boards was studied. 

The results provide evidence that insulation boards and hardboards of good quality 
can be produced from papyrus by relatively simple and inexpensive procedures. 


INTRODUCTION 


The papyrus of Israel grows in the Huleh Valley of the Galilee, a swampy area of 
about 1000 hectares around the Upper Jordan. A quantitative survey carried out 
recently! showed the density and rate of growth of the papyrus to be approximately 
45 tons of o.d. papyrus per hectare per annum. The soil on which the papyrus grows 
is wholly composed of deep deposits of peat, and only at a few isolated spots, 
especially on the banks of the Jordan, has some sandy soil accumulated and a vege- 
tation of papyrus mixed with Phragmites communis and Polygonum developed. 

The papyrus plant is composed of three main parts, differing in their physical 
structure and chemical composition: the efflorescence (or head of the stalk—see 
Plate 1), the rind and the pith. As the averages of 50 determinations on single stalks 
collected in the Huleh area, the following figures were obtained in this laboratory 


* This paper was originally prepared as a contribution to the International Consultation on Insula- 
tion Board, Hardboard and Particle Board, convened under the auspices of FAO and ECE in Geneva, 


January 21—February 4, 1957. 
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for the percentages of these three parts of the oven dry weight of the whole papyrus 
plant: efflorescence — 33%; rind — 37.5% and pith — 29.5%. It seems that the 
Israeli papyrus is markedly different from that of Northern Rhodesia, which is 
said to contain 80% of pith3. There seems also to be a considerable difference in the 
fibre structure between Israeli papyrus and that of Lake Tana and Gabon. While 
the mean fibre length and width of the Israeli papyrus was found to be 0.87 mm 
and 8.65 yw respectively, those of the Tana and Gabon papyrus were reported to be 
1.75 mm and 12 p45. 

Although papyrus seems to have been one of the first raw materials from which 
a kind of paper was made, and was widely used in ancient Egypt and to a consi- 
derable extent in the Hellenistic and Roman empires, its use was discontinued in 
the Middle Ages, and only in the Twentieth Century has it attracted some interest. 
Several papers have been published on its chemistry and possible utilisation. A 
partial review of these investigations is included in the introduction to the paper | 
of Wise and El-Taraboulsi3. The alkaline cooking of papyrus has been tried by Levi 
and Marini+, with a yield of 39% of unbleached and 35% of bleached low quality 
pulp obtained on cooking with 18% alkali at 155°C for four hours. Badawi tried 
a mild cooking procedure, using 10-15% alkali for various times and temperatures, 
obtaining slow pulps with strength properties comparable to those of straw. The: 
soda pulping of papyrus was also tried by Beam’, but the possibility of producing 
insulation and hardboard from papyrus seems not to have been sufficiently investi- 
gated. Only in Uganda is the papyrus utilised to a certain extent for the production 
of cardboard8, | 

In Israel the papyrus is used only to a small extent for the manufacture of mats, 
although some experiments were carried out during World War II on the production 
of insulating panels for cold storage, using partly peeled papyrus chips impregnated 
with wood tar9. | 

Since it has been found that the quantities of papyrus in Israel are not sufficient 
for the production of pulp, work was undertaken in these laboratories on the possible’ 
utilisation of the papyrus plant for the production of insulation boards and hard- 
boards, so that it could serve as a part of the furnish of the Sefen Ltd. company in 
the Jordan Valley. As preliminary experiments showed that the presence of the pith 
does not present any serious difficulties in the production of the wallboards, the: 
investigation was carried out on the whole and undepithed papyrus stalks, but with 
the exclusion of the roots. The chemical composition of the various parts of the 
papyrus stalks and the full data on the dimensions of the fibres, determined in con- 
nection with this work, will be published elsewhere. | 


EXPERIMENTAL PROCEDURE 
The papyrus used in the experiments was brought to the laboratory from the Huleh 
Valley stands and chopped in a hay-cutter into pieces of approximately 4-3” length. 
The moisture content of the air dry material (after several days’ storage in the labo- 
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ratory) was 17%, while the green material contained over 75°% moisture. Quantities 
of 1200 or 900 grams oven dry papyrus were used in the single experiments. 


A rotary digester of 14 liters capacity, electrically heated and equipped with a 
pressure gauge, thermometer and thermostat, was used in the experiments. Since the 
construction of the digester did not allow steam to be fed during rotation, the 
cooking was carried out in the following way: 

The digester was electrically preheated to 120°C, the papyrus and the cooking 
liquor were introduced, and direct steam applied for 7-8 minutes until the desired 
cooking temperature was achieved; the steam was turned off and the electrical 
heating and the rotation of the digester were resumed. The cooking time at pressure 
was usually 20 minutes, and during this time direct steam was introduced 2-3 times 

for periods of 2-3 minutes. At the end of the cooking the pressure was released and 
the contents removed from the digester and weighed in order to determine the amount 
of the condensed steam. The pH of the cold water extract of the cooked chips was 
usually determined with a glass electrode. Samples of the cooked chips were analysed 
after the adhering liquor was removed by a thorough washing with cold water. 


The cooked papyrus was refined at 3000 revs/min in the 12” laboratory single 
rotating disc refiner designed in these laboratories (see Plate 3). The freeness of the 
refined stock was determined in the Williams precision freeness tester on a suspension 
of 4 grams. The stock was collected in a stationary wire box. For the determination 
of yield the stock was centrifuged up to 80% moisture, weighed and the moisture 

contents determined on several samples. 


While for the preparation of the hardboards the refined stock was screened on a 50 
mesh screen for partial removal of the fines, the insulation boards were made from 
the unscreened material. 

A quantity of refined stock correponding to 350 grams of oven dry material was 
taken for the preparation of a single board with the dimensions of 30x 301.2 cm. 
The boards were prepared on a sheetmaker with a wire of 25 mesh. The wet boards 
were cold pressed in a hydraulic press between two wire nets at predetermined pres- 
‘sures. The cold-pressed boards were dried at 150°C in a drying oven with cir- 
culating hot air. The insulation boards were conditioned for 24 hours in a constant 
‘temperature (20°C) and relative humidity (65%) room and thereafter tested. 

For the preparation of the hardboards the cold pressing (see Plate 4) was carried 

out at 100 kg/cm? pressure, so that the moisture content was reduced to 50%. 
The hardboards were formed by pressing at a temperature of 200°C, measured 
with a copper-constantan thermocouple on the surface of the board, and at 100 kg/cm? 
pressure. The time of closing and of reaching this final pressure was 30 seconds. 
In order to avoid ‘“‘water burns” and blisters a breathing period of 30 seconds at zero 
pressure was introduced. After reclosing, the final pressure was reached in 15 seconds 
and the hot pressing was continued for 6 minutes. The nominal thickness of the hard- 
boards was 1/8”. In some of the experiments 0.6% rosin size and 1% alum were 
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added to the refined stock before the sheet formation. In others, up to 4% of linseed 
oil was added to the stock after fixing the rosin size with the alum. 


Some of the boards were tempered by heating in circulating hot air at 180°C for 
2 hours or by brushing 3% linseed oil on the surface of the boards while hot 
from the press and heating similarly for 2 hours. 


A full-scale mill trial was carried out at Sefen Ltd. in the Jordan Valley, a combi- 
nation insulation and hardboard plant with equipment supplied by the U.S. 
Wallboard Machinery Co. Four and two-tenths tons of papyrus stalks containing 25% | 
moisture harvested in the Huleh Valley, were chopped into 1-2” length pieces and fed | 
into a 65 m3 spherical digester heated with direct steam. During the feeding, the 
papyrus chips were sprayed with 75 kg sodium hydroxide dissolved in 2100 litres 
of water, so as to have a liquor ratio of 1:1. The cooking conditions were: raising the | 
temperature to 130°C in 11 minutes; additional raising to 150°C in 4 minutes; time at 
150°C — 20 minutes; release of pressure — 10 minutes. After dumping, the cooked 
papyrus chips were refined on Bauer refiners. The refined pulp was screened on a 
vibrating screen of 1.5 mm slots, and the sheet formed at a speed of 10 feet per | 
minute on an Oliver forming cylinder. The board mat before entering the Coe | 
dryer had a moisture of 70%. The drying of the insulation board was carried out | 
at a temperature of 165—170°C. 


TESTING METHODS 


The insulation boards were tested according to the British Standard! specification. 
The modulus of rupture was determined with the Defibrator flexural resistance | 
tester on board specimens of 100 mm span and width. On several insulation boards | 
the deflection on break was measured during this determination. | 
The tensile strength was determined on the Lorentzen and Wettress tester, using» 
board strips of 100 mm width and span. The results are given in kg/cm2. 

Onseveral samples the brightness was measured by means of a Photovolt instrument. 

The specific gravity was calculated from the dimensions and oven dry weight of 
the boards. 

In place of the standard apparatus for the heat conductivity measurements, an 
apparatus was used in which the temperature gradient was 80°C (from 100°C to 20°C). | 
A dynamic method of measurement was adopted. The temperature of the hot side | 
was kept constant at 100°C by a copper vessel containing boiling water and the| 
bottom of which was in contact with the board, and the rise in temperature of an | 
insulated copper cylinder in contact with the other side of the board was measured | 
by means of a differential thermocouple and a sensitive galvanometer. The heat 
conductivity was calculated from the change of the temperature difference with 
time, according to the following equation: | 


_ log(D./D) 2.303 M-S-h | 


Ro ee ee 
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where: Do = Deflection of the galvanometer, proportional to the temperature 
difference at zero time. 
D = Deflection of the galvanometer at time ¢ minutes. 
K = Heat conductivity of the board. 
M = Weight of heat sink (copper cylinder) = 370 grams. 
S = Specific heat of copper = 0.093 cal/gr. 
L = Thickness of board specimen in cm. 
A = Area of contact = 15.2 cm2. 
The deflection D was measured at time intervals of 3 minutes and K was calculated 
as an average of 6 measurements. 
Since the heat conductivity depends on the temperature, it is realised that the 
above method yields only relative values. 
The surface hardness of the hardboards was measured according to the procedure 
of Brinelll!. A steel bearing ball of 10 mm diameter was pressed on the hardboard 
| in the Hinde and Daush crush tester at a load of 100 kg for 15 seconds, and the 
diameter of indentation was measured 60 seconds after releasing the pressure. The 
surface hardness was expressed in kg/mm? using a table given by Kollman!!. 
_ The puncture resistance of the hardboards was measured in the Beach tester!2. 
The range of this tester was found suitable for the boards prepared in the present 
| investigation. 
_ The water absorption of the hardboards was determined by immersing 100 x 100 
specimens under 1”’ of water for 24 hours. 


PREPARATION OF INSULATION BOARDS FROM PAPYRUS 


( Preliminary experiments on the laboratory production of softboards from whole 
papyrus stalks are summarised in Table I. The liquor ratio in these experiments 
1 was 1:0.92 at the beginning of the cooking and 1:1.78 at the end. In the 
' cases in which 900 grams of papyrus containing 17% moisture were taken, 
| 540 ml of water in which the alkali or the bisulphite were dissolved were added to 
the digester. The liquor ratio changed owing to the condensation of the direct steam 
| during cooking. It is seen from Table I that the specific gravity and the tensile strength 


of the boards increase with an increase in temperature. The final pH of the liquor 
| 


was high and decreased with an increase in temperature, indicating a more efficient 
| utilisation of the alkali. When the temperature was increased to 150°C the final pH 


dropped to 7.5-8 (see Table II). In the cooks with sodium bisulphite and in those in 
which no chemicals were added, the pH was around 6, suggesting that during the 
i digestion of the papyrus the drop in pH is caused by the consumption of alkali and 
) not by the formation of organic acids in appreciable amounts. In the case of similar 
| wood digestion the final pH of the spent liquor is much lower !3. 

| The yields of the pulp obtained in the various cooks seem to be approximately 80%, 
| which is significantly lower than the yield obtained from hardwoods and tropical 
| woods!3, This lower yield might be due to the relatively high amount of water soluble 
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matter in the papyrus (12-13%), as compared with wood. The results of Table II 
indicate that the yield is not noticeably influenced by the chemicals used for the 
digestion and that a steam and water treatment produces the same yield as alkali 
and bisulphite cooks. Since the specific gravity and the strength properties of the 
boards obtained on using alkali or bisulphite are not significantly different from 
those obtained when cooking with water and steam only (see Table IID), it is deduced 
that insulation board can be advantageously produced from papyrus on cooking 
at 150°C without the addition of chemicals. 

In Table IV the influence of the liquor ratio in the digesting stage on the properties 
of the insulating boards is summarised. It seems that the amount of water added to 
the digester is critical, and that below an initial ratio of 1:0.92 the strength is unsatis- 
factory. In the “‘dry’’ steam digestions (exp. 260 and 270), the material had not been 
sufficiently softened and the pulp obtained was hard and very coarse, and a higher 
energy was required for refining it. The pulp of experiment No. 280 was softer and more 
easily refined, but the strength was still low, while in the experiments with a higher 


liquor ratio the pulp and the boards obtained were satistactory. The liquor ratio | 


of 1:0.92 was therefore adopted in all laboratory experiments reported in this paper. 
Several series of experirnents were carried out in order to assess the dependence 
of the properties of the insulating boards on the mechanical variables. It is seen from 


Table V that while the strength and specific gravity of the insulation board both | 


increase with the increase in the number of times it is passed through the refiner, 


the most significant change in properties occurred between the boards which were | 
refined with one or with two passages through it. The distance between the discs was | 


identical in all the experiments, and the data of Table V therefore give a rough 
estimate of the influence of the degree of refining, showing that two passages of that 
material through the refiner are necessary for the preparation of a satisfactory in- | 
sulating board from papyrus. 

The pressure used in the wet pressing stage seems also to influence strongly the 
properties of the insulation boards. It seems that the use of pressures in the range 
of 5-50 kg/cm? yields boards having properties within the specifications required 
(Table VI). 


The properties of the laboratory-made insulation boards are compared in Table VII | 
with the properties of the boards obtained in the mill trial described above (see) 
Experimental). The similar characteristics of both boards suggest that the close) 
setting of the Bauer refiners used conformed closely to the two-passage refining | 
procedure employed in the laboratory. A coarser refining in the mill resulted in) 
boards having a considerably lower strength. The commercial board tested was taken | 
from one of the batches produced from a furnish consisting of 75° unbarked| 
Eucalyptus rostrata chips and 25% ‘okume’ wood waste at Sefen Ltd. by essentially) 
the same procedure as used for the papyrus, except that the cooking temperature} 
was 180°C. The papyrus insulation board seems to be stronger and brighter than\ 
the eucalyptus board and the thermal conductivity considerably lower. The thermal | 
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conductivity of the eucalyptus board when measured in the standard apparatus 
according to the British specification was 0.46 BTU/h/F2/°F/inch!4. The heat 
conductivity of the papyrus board, being lower than that of the eucalyptus board, 
seems, therefore, to be well below the upper limit of that property as demanded by 
the British specification. 

No marked difficulties due to the presence of pith were encountered in the present 
investigation. The main objections to the presence of the pith arising in the production 
of pulp and paper from pith-containing raw materials hardly arise in the case of board. 
Since no chemicals or very small amounts of them are needed in the cooking stage 
of the board, according to the present investigation, no increased consumption of 
chemicals will be caused owing to the greater surface area and ease of penetration 
through the pith. The higher absorbency of the pith also does not play any significant 
part in the various stages of the board production owing to the rather limited require- 
ments for purity of the pulp. On the other hand, the freeness of the pith-containing 
papyrus pulp was found not to be seriously affected, as the speed of running on the 
board machine and the rate of draining were satisfactory, and not significantly 
different from those of wood pulps used in the same conditions. 


An interesting feature of the papyrus insulation boards is their deflection at 
break, which is significantly higher than in the case of the eucalyptus boards. It 
seems to be closely linked with the higher strength of the papyrus wet fibre mat 
observed when descending from the Oliver forming machine. It is suggested that 
the difference in the fibre length distribution in the raw material and in the size of 
the fibre bundles of the refined stocks is responsible for these differences. This sug- 
gestion is substantiated by the screen analyses summarised in Table VIII. It is seen 
that the papyrus board consists generally of longer fibre bundles and it contains a 
considerably lower percentage of material passing through the 100-mesh screen. 


In Figure 1 the fibre length distribution of the pith-containing papyrus stalks is 
compared with that of a representative sample of the E. rostrata as used in the Sefen 
Works Ltd. The length measurements were carried out on 400 fibres in each case. 
The material was prepared by macerating the wood or papyrus chips with a mixture 
of equal volume parts of glacial acetic acid and 100 vol. hydrogen peroxide and heating 
in a water bath until disintegration on vigorous shaking. The percentages given 
in the figure are percentages of the cumulative length of all fibres measured. It is 
seen from Figure 1 that while the longest eucalyptus fibres are in the range of 
1.1-1.2 mm, in the case of papyrus more than 9% of the weight of the material 
consists of fibres in the range of 1.2-1.9 mm. Although no definite data on the rela- 
tionship between the fibre length distribution of a given raw material and the fibre 
classification of the refined stock prepared from this raw material are available, 
it may be assumed at least qualitatively that the latter is largely influenced by the 
former. It may thus be expected that a papyrus of the kind measured by Vidal>, 
having fibres of 1.7 mm average length, might yield boards of a still greater flexibility. 
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PREPARATION OF HARDBOARDS FROM PAPYRUS 


The pulp for the preparation of the hardboards was prepared in exactly the same 
way as described above for the insulation boards, except that before the formation 
of the boards the pulp was washed on the Somerville classifier over a 50-mesh screen. 
The removal of fines amounting to 5% calculated on the oven dry papyrus stalks 
by this screening treatment was necessary, since the boards which were prepared 
from pulp without screening developed burns and water pockets during pressing, 
in spite of the breathing period used. The freeness of the pulp after this screening 
was 27 Williams seconds (620 ml Canadian Standard Freeness), while without 
screening the freeness was 34 seconds (575 ml). 


The conditions described above (see Experimental) for the pressing of the hard- 
boards were the optimum conditions chosen after a long series of trials. Pressing 
at higher temperatures than 200°C caused ignition of the boards in the press, whereas 
lowering the temperature and pressure yielded weaker boards. 


In Table IX the properties of some of the boards prepared from the unsized stock 
are summarised. It is evident that except for water absorption the properties of these 
boards are satisfactory and close to the British Specification for standard hardboard 
of 1/8’ thickness. Heat tempering does not change significantly the strength of the 
board, although it lowers the rather high water absorption. The oil tempering followed 
by heat treatment markedly increased the strength, but did not improve the water 
absorption beyond that of heat tempering alone. Increasing the amount of oil did 
not yield better results. 

Sizing of the pulp with rosin and alum appreciably increased the modulus of 
rupture of the board (see Table X), but only a tempering treatment with linseed oil 
reduced the water absorption of these boards to 25%, thus placing them within 
the British Specification for standard hardboard. 


It is noteworthy that the inclusion of 4 & linseed oil in the sizing treatment increased 
the water absorption of the untempered boards (compare boards 340 and 340b in 
Table X). Heat tempering of this board, however, reduced the water absorption 
considerably. 

The surface hardness as measured by the Brinell method was not noticeably 
influenced by the treatments described in Table X, and was in all cases 7.6 kg/mm2 
(see Table XI), indicating that the papyrus board is fairly hard 11, 


The puncture resistance of the papyrus hardboards as measured in the Beach 
tester was similar to that of the commercial boards (see Table XII). when considering 
the differences in thickness. 

An overall comparison of the properties of the papyrus hardboards described in 
the present paper with those of several commercial hardboards tested in these 
laboratories is presented in Table XI. It seems that the quality of the papyrus hard- 
boards is very similar to the quality of the Dutch hardboard produced from straw, 
although the water absorption of the papyrus boards is lower and approaches that 
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of the Swedish-made wood hardboard. The colour of the papyrus hardboards is 
gray, and its brightness is close to that of the wood hardboard. The papyrus hard- 
boards can thus be classed as standard hardboard, while Insulite seems to be 
a super-hardboard and the “commercial building board” is a medium hardboard. 


The difficulty encountered in this investigation in obtaining super-hardboard 
from the papyrus may be ascribed to its low lignin content. It is possible that an 
addition of a small amount of a suitable resin will considerably increase the strength 
properties. Experiments in which the papyrus pulp was blended with varying amounts 
of hardboard-grade pulps made from deciduous woods yielded super-hardboard. 


SUMMARY AND CONCLUSIONS 


1. The production of insulation board and hardboard from undepithed Israeli 
papyrus was investigated under laboratory conditions. 


2. When cooking the papyrus with 3% alkali at a 1:0.92 liquor ratio, the final pH 
decreased and the specific gravity and tensile strength of the boards obtained increased 
with an increase in temperature. 


3. Good quality insulation boards can be obtained when cooking with water at 
150°C, and no advantage was found in adding up to 2% alkali or sodium bisulphite 
of the weight of the o.d. papyrus. 


4. The yields obtained were approximately 80% 0.d. The pH in the cooks with 
water was 6.0, while in the cooks with alkali it ranged from 7.5 to 10.6. 


5. The liquor ratio in the cooking stage is critical for the strength properties of 
the insulation board, and only at ratios of approximately 1:1 can boards of satisfactory 
properties be obtained. 

6. The degree of refining and wet pressing strongly influences the specific gravity 
and strength of the insulation boards. 

7. The heat conductivity of the papyrus insulation boards is satisfactory and is 
not appreciably influenced by cooking conditions. 

8. The flexibility and wet strength of the fibre mat appear to be significantly higher 
in the case of papyrus than in the case of eucalyptus. This difference seems to be 
related to the differences in the fibre length distribution of the respective raw materials 
and to the screen analysis of the insulation boards. 

9. Ina full-scale mill trial, insulation boards from papyrus were obtained having 
properties within the range of the British Standard specifications, similar to laboratory 
samples. 

10. Hardboards were produced from papyrus by hot pressing of wet board mats at 


a temperature of 200°C and a pressure of 100 kg/cm? for 6 minutes after a breathing 


period of 30 seconds. 
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11. Screening of the hardboard pulp through a 50-mesh screen resulted in a loss of 
approximately 5% fines and decreased greatly the development of burns and water 
pockets during the pressing stage. 


12. Heat and oil tempering lowered the high water absorption of the papyrus 
hardboards but only oil tempering increased the strength of the board. Sizing with 
rosin and alum followed by oil tempering yielded hardboards with a water absorption 
oF 255%: 


13. The papyrus hardboards conformed to the British specification for standard 
hardboard and had properties similar to those of commercial Dutch and Swedish 
hardboards. 


14. It is proved by the present investigation that insulation boards and hardboards 
of good quality can be produced from papyrus by relatively simple and inexpensive 
processes. 
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TABLE I 
Pulp and board from papyrus: preliminary experiments 
Initial liquor ratio = 1:0.92; final liquor ratio = 1:1.78; NaOH = 3% on o. d. papyrus; time | 
schedule = 8 + 20 min; one passage through refiner 


No. of Grams of o. d. Temperature Final Sp. gr. of Tensile 


exp. papyrus taken of cook °C pH board strength 
per cook gr/cm3 kg/cm2 
110 1200 110 0.26 4.8 
120 1200 130 10.6 0.29 6.5 
130 900 140 9.7 0.32 16.5 
140 900 140 9:2 0.35 20.3 
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TABLE Il 


191 


Time schedule = 8 + 20 min.; Temp. = 150°C; Initial liquor ratio = 1:0.92 


No. of Chemical &% on o. d. Yield o. d. on o. d. Final 
exp. papyrus és pH 
210 2.0 — NaOH 80 8.3 
220 1.5 — NaOH 74 7.6 
230 2.0 — NaOH 79 Tes 
240 Water only 80 6.0 
250 2.0 — NaHSO; 79.5 5.9 


TABLE Il 


Specific gravity and tensile strength of papyrus insulation boards prepared from pulps cooked with 
different chemicals 
Temp. = 150°C; initial liquor ratio = 1:0.92; pressure used in wet pressing = 50 kg/cm2; one 
passage through refiner 


No. of 
exp. 


220a 
220b 
240a 
250a 


Chemical % on o. d. 


papyrus 


1.5 — NaOH 


1.5 — NaOH 
Water only 


2.0 — NaHSO, 


TABLE IV 


Tensile strength 


kg/cm? 


Influence of the water—to-—papyrus ratio in the digester on the properties of the insulation boards 


Time schedule = 8 + 20 min.; temp. = 150°C; two passages through refiner; pressure used in 
wet pressing = 50 kg/cm2; amount of papyrus per cook = 750 gr o. d. or 900 gr a. d. 


No. of Water added 
exp. grams 
260 0 
270 250 
280 400 
240b 540 
290 750 


Intial ratio o. d. 


papyrus/water 


1:0.2 
1:0.4 
1:0.74 
1:0.92 
eZ 


Sp. gr. Tensile 
gr/cm3 strength 
0.250 2.4 
0.260 5.0 
0.250 4.3 
0.280 15.4 
0.290 10.8 


TABLE V 


Properties of insulation boards prepared from papyrus with varying numbers of passages through the 


refiner 


Cooking schedule = 8 + 20 min.; temp. = 150°C; pressure used in wet pressing = 50 kg/cm?; 
initial liquor ratio = 1:0.92; normal thickness of board = 12 mm 


No. of passages 
through refiner 


1 
2 
3 


4 


Sp. gr. 
gr/cm3 


0.26 
0.28 
0.29 
0.32 


Tensile strength 
kg/cm2 


Modulus of rupture 
kg/cm2 


23 
31.6 
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TABLE VI 
Properties of insulation boards prepared from papyrus at various pressures in the wet pressing stage 
Cooking conditions as in Table V; two passages through refiner 


Pressure Sp. gr. Tensile strength Modulus of rupture 
kg/cm? gr/cm3 kg/cm? kg/cm? 
1 0.23 3.0 
5 0.27 14.0 20.2 
50 0.28 15.4 23.0 
150 0.37 33.0 36.1 
TABLE VII 


Properties of the insulation boards made from papyrus in the mill scale trial, as compared with laboratory 
boards and commercial Eucalyptus rostrata boards 


Relative ther- 
Kind of board Thickness Sp. gr. Deflection Modulus of Brightness mal conduct- 


at break rupture ivity BTU/h/ 
mm gr/cm3 mm kg/cm2 %G.E. F?/F/inch _ 
230b (lab.; 2% 
NaOH) 10 0.29 14 19 19 0.86 
240b (lab.; water) 11 0.28 14 23 23 0.82 
Mill trial 
(fine refining) 12.4 0.28 10 23 19 0.82 
(coarse refining) Veoh 0.28 11 19 
Eucalyptus 
comm. board 12.6 0.30 7 132 12 1.03 
TABLE VIM 
Fibre classification of papyrus and commercial boards (in per cent of dry weight) 
Kind of board Above 12—18 18—30 30—S0 50—100 Below 
12 mesh mesh mesh mesh mesh 100 mesh 
Papyrus (mill trial) 10.1 16.8 11.9 19,3 22.0 19.9 
Comm.eucalyptusboard 8.1 9.0 123 17.2 24.0 29.4 
TABLE IX 


Properties of papyrus hardboards prepared from unsized stock with and without tempering 


Pulp: as for insulation board; cold pressing = 100 kg/cm2; temp. of hot pressing = 200°C; time | 


cycle = 30 sec. to pressure; 30 sec. breathing; 6 mins. at 100 kg/cm? pressure 


No. of board Thickness Sp. gr. Modulus of rup- Water absorption 
mm. gr/cm3 ture kg/cm2 Vi, 
320; untempered Qi 1.03 338 95 
320a; heat tempered 
(180°C, 2h) 2 1.07 337 67 
330; untempered Ses} 1.05 320 82 


330a; oil tempered; 

brushing with 3% 

linseed oil + heat 2.7 1.07 360 68 
330b; oil tempered 

(4% linseed oil) 3h} 1.07 350 67 
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TABLE X 


Properties of papyrus hardboard prepared from sized stock with and without tempering 
Conditions as in Table [X 


No. of Sizing treatment Tempering Thickness Sp. gr. Modulus of Water 
Board gr/cm3 rupture absorption 
mm kg/cm2 Ye 

330 Nil Nil 33 1.05 320 82 
340 0.6% rosin+1% Nil 

alum 35) 1.06 370 64 
340a 0.6% rosin+ 1% Heat 

alum 3).5) 1.07 390 MS 
340b 0.6% rosin + 1% 

alum + 4% linseed oil Nil 3.4 1.07 360 85 
340c 0.6% rosin+1% 

alum + 4% linseed oil Heat 3.4 1.07 390 38 

TABLE XI 


Comparison of laboratory produced papyrus hardtoards with commercial wood and straw hardboards 


Kind of hardboard Lab. papyrus Straw, Dutch! Wood, Insulite, duo- Commercial 
Nivalite Swedish! faced, Finnish building board, 
Softwood 
Thickness, mm 
nominal 33.5) 3.86 325) 3.14 4 
measured S55) 4.1 3.40 3.4 4 
Specific gr., gr/cm3 1.07 1.06 0.96 1.20 0.96 
Modulus of rupture, 
kg/cm2 based on 
nominal thickness 398 335 650 234 
measured thickness 390 345 305 580 234 
Water absorption 
% 24h 25 58 23 22 35 
Brinell surface 
hardness, kg/mm2 7.6 7.6 Se 10.7 5.9 
Beach puncture re- 
sistance, in. OZ. 570 690 620 750 
Brightness (green 
filter) 7, GE. 13 24 13 14 25 


1 The samples were kindly supplied by Defibrator AB, Stockholm. 
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Figure 1 


Fiber length distribu- | 
tion of papyrus and | 


ee 53 ae Eucalyptus rostata 


Figure 2 
Papyrus in the Huleh Valley; note the efflorescence. 


Figure 3 
Papyrus, general view of fibrous mate- 
rial; note fibers, characteristic vessels 
and parenchyma cells. x 220. 


Figure 4 
The single rotating disc refiner. 
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Figure 5 | 
The hydraulic press, with thermostatically controlled and electrically heated plates. | 


ON THE MEMBRANE STATE OF STRESS IN A CURVED CIRCULAR TUBE 


A. KORNECKI 


Department of Aeronautical Engineering, Technion—Israel Institute of Technology, 
Haifa 


The computation of stresses in a thin-walled shell is in general a rather complex 
problem and to date exact solutions (within the accuracy limits of Love’s theory of 
shells) have been obtained in only a few cases. 


For most engineering purposes, however, it is permissible to apply the solution 
based on the membrane theory, superposing if necessary the stresses of the “‘edge 
effect”? on the membrane stresses in order to satisfy the boundary conditions. 


The toroidal shell, as is well known, has the singularity that computations based 
on the membrane theory result in indefinite values of stresses and/or strains at points 
on the lines 9 =0, 2, along which the average curvature of the shell’s midsurface 
changes its sign. Since parts of toroidal shells are encountered quite often in practice, 
it is of interest to ascertain whether a membrane state of stress is at all possible 
in a shell of this type under any load. 


Consider for this purpose a circular curved tube (Figure 1) loaded at its end 
section y = 0 by arbitrary unit forces Tp(@) and S(Q) and supported atp = Y » Say 
(without loss of generality) by means of a diaphragm, rigid in its own plane and 
providing no restraints against displacements perpendicular to its plane. 


According to the membrane theory, we may write the equilibrium equations of the 
toroidal element under investigation in the following form: 


1 le cos @ 1 as 
fo Rash) O81, RGA Aan 6) - ep | 
1 os cos @ 1 Ty 
ee ee ee =O, 1 
(0 te ROL) ORT ein 6). Op (1) 
To Ty sin @ | 


gor aR aay a 
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Figure 1 
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Some of the symbols used in these equations are indicated in Figure 1, the others 
being self-explanatory. The parameter / equals 


r 


Expanding the left parts of Eq. (1) in power series with respect to the small parameter 
A we obtain 


Tg =—A(I —Asin @ + /?sin2O —...)Tysin@, 
To as ; 
36 = 41 (Te — Te) cos © — ace! (1 —Asin O + A’sin?@ — Ory (3) 
os Ty 
se = — 42S cos 6 + ae (1 = Asin Ost A SiO een ye 


Let us seek the unknown stress resultants in the form of power series, too: 


T, (8, , 4) = Ty (0,9) +4Ty (0,9) +#T2(O.9) ++--1= O03 py 
S(O, 9, a) = So(O,~) + 4S; (0,9) + 4752 (0,9) +--+ 


where 7;), T;,,..- So, S,...do not depend upon A. 
Substituting series (4) in Eq. (3) and comparing terms with equal powers of A we find 


Te 00) G =O 7 = 01,0...) ’ 

Bat 2597 sino 29 Pe sino TR (S21 a aia OF: ‘i (5) 
or 

T= 0p = 0 Sac (asin On. .2 (6) 


where C is a constant of integration. 

Expressions (5) have been obtained by consistently setting equal to zero terms 
which are indefinite at 9 = 0 or/and 0 =z. 

Making use of the relations known from the theory of shells we state that to 
stresses (6) correspond, for instance, the following displacements 


2Cr(i + ; d@ 
Ug = 9, ie == (0), up = ST (+ Asin @) SOF isin ©’ (7) 


which are definite at all © and fortunately satisfy the boundary conditions at p= 9,. 

Accordingly we conclude that if the unknown functions 7;, S admit of expansion 
in power series (4) the only possible state of stress is defined by (6). Hence 
we may load the end section p = 0 by forces JT, S defined by (6) and statically 
equivalent to a twist couple, as may be readily checked. Any other edge forces 
cause bending moments, at least in the regions adjacent to the singular lines 
0 =0,2. 


200 A. KORNECKI Bull. Res. Counc. of Israel 


Similar results were obtained by Sh. S. Mekhovrishvili!, seeking the stress 
function W corresponding to the full torus (y, = 2) in the form of infinite series, 
W= ZX OG (Oe, (8) 
k=-a 


by a much more involved mathematical analysis. 
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LETTERS TO THE EDITOR 


‘Note on the paper “Ejector Refrigeration from Low Temperature Energy Sources’!, M. VAN MENTS, 
‘The Research Council of Israel, Jerusalem. 


Mr. H. Tabor has drawn my attention to the fact that by raising the boiler pressure Po, the efficiency 
of the cycle decreases if use is made of the results of eqs. (2) and (3) of the above article}. 


It is evident, however, that the authors have misunderstood a very important point in the article 
of Work and Haedrich? to which they referred. The value Po/P-=1/7.3 (p. 467)? refers to the 
minimal suction pressure (at zero entrainment) attainable and this occurs only if Ps = 5.5 Pe. 


Consequently it is permissible to use eq. (2) of Mizrahi et al.! in this special optimum case 
only and then (P-/P»)°-85 = constant = 0.24! 


i 
r 
tf 


As a result, all that follows after eq. (2) of Mizrahi et al.! is quite erroneous. One could make 
/ an attempt, however, to use the results of Work and Haedrich? by applying the following procedure: 


From Figure 4 in their article, it follows that 


Po 1) Pel— + Po, (1) 


If this relation held for Freon 12 (although this is by no means certain)*, then 


on condition that Ps > Py > P-/7.3. 
Po © 162 — Po (2) 
! 


| in the case where T- = 86°F (Pe = 108; Py > 15). 


_ By inserting the above eq. (2) into eq. (1) of Mizrahi et al.!, one obtains (if Py = 27) 


(3) 


0.313 i pee Ue ( >). 


tue ky 
Wo Py — 220 


on condition that 540 < Ps < 600. 
The condition shows that the term (108/Ps)°:85 of the above eq. (3) hardly depends on Ps, so 
for this eq. (3), one could also write 

— 540 


elena 4 
big 7m — 220’ @) 


also on condition that 540 < Ps < 600. 


* Other constants will also do, and the approach remains the same. 
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It can immediately and very clearly be seen from the above eq. (4) how the cycle efficiency 
increases with the boiler pressure. 
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Photoelastic analysis of plates by bonding polaroid foils between plastics, Cu. H. LERCHENTHAL, Division 
of Structural Engineering, Technion—Israel Institute of Technology, Haifa. 


A number of methods are known for the photoelastic analysis of plates in transverse bending!. 
I have now cemented polarising foils (a) (see Figure 1) between two sheets of the same photoelastic 
material and of equal thickness (b)?. The advantage of such a procedure is the possibility of a separate 
analysis of both ‘‘half’-plates (b), thus permitting the determination of bending and membrane 
stresses. An ordinary transmission polariscope can be adapted to this method by removing alternately 
polarizer or analyzer. Preliminary tests indicate that this method is feasible. 


Natural 
Light 


Figure 1 


The polarizing foils are sheathed in cellulose acetate butyrate. Highly stress-sensitive plastics 
like Catalin Resin 61/893, Columbia Resin CR-39, Araldite, Castolite, or WP-1527, which are cross- 
linked polymers of the phenol-formaldehyde, allymer, ethoxylene, or polyester families respectively, 
are difficult to cement to the polarising foils. Upon the advice of Dr. Z. Rigbi of the Rubber Research 
Association, Vulcafor VCC (a polyisocyanate produced by ICI Ltd., London) was tried with success, 
giving ultimate strength values of up to 210 kg/cm? in shear. Araldite Cement 103 with Hardener 951 
(products of CIBA A. G., Basle) was found to be somewhat weaker (shear strength values of up to 
85 kg/cm?), although considerably easier to manipulate. 
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/ Some thermodynamic properties of Freon-12, D. PNUELI, Technion — Israel Institute of Technology, 
Haifa. 


The customary equation of isentropic change of state, Pv“ = const. (where x = Cr/cev) is ultimately 
based on the assumption of an ideal gas; it is necessary to introduce an exponent of the isotherm x7! 
i defined by Pye = const., where 


Vv 
~ PO/op), (1) 


0 


{which is unity for an ideal gas. The isentropic exponent would then be written x. and be determined 
iby 


Hs = XT Col v. 


ur and xs are, by definition, functions of any pair of thermodynamic properties and can be repre- 
jsented accordingly, as shown below. 

| With the publication of the new tables for dichloro-difluoromethane, CC1,F, (Freon-12)2, it 
» became possible to calculate the two exponents for this substance. An attempt to follow an earlier 
‘publication for Freon—223, by establishing an empirical equation for cy = f(P;T), from which the 
; Clausius equation would yield v = f(P;T), failed, due mainly to the very wide range of data included 
in the tables, up to slightly above the critical pressure — nearly three times the range of the tables 
‘for Freon-22. 

It was necessary, therefore, to proceed differently, and all the calculations were done on the basis 


ifor steps of 10°F; similarly, (Av/AP)r and (Av/AT)p replaced the differentials for the determina- 
tion of xr (eq. 1) and of cy, 


(av /eTy., 
~ 1 GaP yn @) 


Cp — Cy = 
The calculations were carried out with an 11”’ slide rule. Intermediate results were submitted to Du 
Pont for consultation, and graphs for cp, cv and cy/cv obtained from the latter showed agreement 
)within less than 1%. Thus, while the method of finite differences may not always be reliable, it 
does not seem to introduce noticeable errors in the present case. 

) The isothermal and isentropic exponents were then calculated on this basis. The results (Figure 1) 
}confirm the general trend demonstrated earlier for Freon—223 , steam1 and ammoniaS. The range 
‘of the curves corresponds to that of the original Freon tables; the crossing of two sg isobarics may 
‘be due to the fact that one of them represents a pressure slightly above the critical, but the tables 
‘do not extend far enough to obtain a complete picture. 


The results permit direct calculation of the velocity of sound: 


a= 4/ 2, Py (3) 


\This equation does not presuppose an ideal gas but only isentropic compression and expansion 
within narrow limits, as caused by a sound wave. As an example, for 200 lb/sq. in and 305°F, the 
‘specific volume may be obtained from the Freon tables as 0.3100 cu. ft/lb and xs = 1.059 from 
| Figure 1, which leads to 


a = V/1.059 x 32.17 x (200 x 144) x 0.3100 = 551 ft/sec. 


{Received June 10, 1958. 
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Figure 1 


The isothermal and isentropic exponents of Freon—12 as functions of temperature for constant 
reduced pressures 7 = P/Pertt. 


7%, P (ib/sq. in) P (kg/cm2) 
0.0234 0.14 0.00985 
8.37 50 Sho) | 
33:5 200 14.05 
67 400 28.10 
100.5 ' 600 42.15 
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Saad Set 


Determination of CaO in hydrated lime, M. RoLEK, Taasiat Even Vasid; J. BARTA, Technion—Israel 
Institute of Technology, Haifa. 


The determination of CaO in hydrated lime is important in production control and in drawing up 
specifications for this material. No direct and rapid method for the determination of free CaO in 
hydrated lime can be found in the literature. The A.S.T.M. method! determines CaO by 
calculation when mechanical (free) moisture, loss on ignition, CO,, SO3, MgO (total), and CaO 
(total) are known. The Russian method? is along the same lines. 

The method proposed here is based on the reaction: 


CaO + H,0 = Ca(OH), + 15,580 cal. 
Received July 23, 1958. 
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iwhich permits determination of the CaO content in the hydrated lime samples by measuring the 
amount of heat evolved during the reaction, which takes place in a suitable calorimeter. 

| Measurements are carried out by means of a simple device, easily obtainable and sufficiently 
iuccurate for the purpose at hand. 

| The calorimeter consists of a Dewar flask of about 400 cm3 capacity and a thermometer allowing 
veadings up to 0.05°C, which is fitted to the flask through a rubber stopper. Mixing is carried out 
by shaking the fiask, with the thermometer and rubber stopper in place. 

} The hydrated lime sample (10-20 g) is introduced with the aid of a wide-stem funneland a small 
(orush. 

| The water equivalent of the calorimeter is determined with CaO (B.D.H. Laboratory Reagent) 
{ind calculated as follows: 


(278 W + H)—QAt 
At 


C= 


H = heat of solution of Ca(OH), to yield a saturated solution in the calorimeter (in this case 21 cal) 
' Q = quantity of water in the calorimeter, g. (in this case 250 g) 
t = corrected temperature increase (a temperature gradient of 1-2°C between the surroundings 
and the calorimeter results in a temperature difference of 0.1°C/hr). 


> 


ii The CaO content of the hydrated lime sample is then determined in the same manner, using the 
} ollowing calorimeter equation: 


CaO (C+ Q)At— HH” -y 100 


(percentage by weight) 278 W, 


“ vhere W, = weight of the sample in g. 
, Several samples of hydrated lime were analysed with results as follows (percentage of sample 
1s received): 


i 


i CaO CaO 
Sample calorimetric method calc. A.S.T.M. method 

1 6.1 8.4 

27.8 32.2 

3 1.5 2.9 


} It may be assumed that the discrepancy between the method proposed and the A.S.T.M. method 
's due to the dead-burnt lime content of the samples, which would appear in the calculated result 
‘put not in that obtained by the calorimetric method. 

') The suggested method seems especially convenient for production control, and due consideration 
should be given to its inclusion in the Israel standard specifications for hydrated lime. 
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BOOK REVIEWS 


PROBLEMS IN WOOD CHEMISTRY. United Nations Food and Agriculture Organisatior 


Proceedings of the Seventh Session of the Technical Panel on Wood Chemistry, held on 8-13 April, 
1956, in Jerusalem and Haifa, in cooperation with the Institute for Fibres and Forest Product: 


Research, Israel Ministry of Agriculture. Weizmann Science Press, Jerusalem, 1957, 136 pp. 


A compendium of lectures by the foremost experts in the field review here the latest developme: 
in several aspects of research in wood chemistry. 

An article on “Radiation Chemistry and Wood” by Prof. H. Mark, of the Polymer Resear 
Institute of the Polytechnic Institute of Brooklyn deals with the new subject of the influence 
nuclear radiation on the polymerisation of monomers, degradation of high polymers and modificati 
of high polymers. The chemistry, morphology and mechanical properties of the compound mid 
lamella of woods are dealt with by Prof. H. W. Giertz of the Technical University of Norw 
Trodheim. An extensive review on the middle lamella of the bast fibres is given by Dr. M. Lew 
of the Israeli Fibres and Forest Products Research Institute, Jerusalem. The morphology, chemi 
composition, biological stability, chemical reactivity and strength properties are dealt with ur 
separate headings; the influence of the middle lamella on the bleaching, swelling, mercerisat? 
moisture absorption, crimp and Young’s modulus of the fibres to which it adheres is also thoroug 
discussed. A summary and a list of 75 references complete what is probably the most exhaus’ 
critical review of the subject today. Prof. K. Freudenberg of the University of Heidelberg discu 
the problem of a possible carbohydrate-lignin bonding in spruce. A fundamental approach to 
problem of adhesives, mainly natural and synthetic high polymers, which are used in wood techr 9). 
is provided by Prof. H. Mark in his excellent article ““A Few Principles of Adhesion,” in which ~ 
theoretical aspects of adhesion are discussed from both the thermodynamic and the molecular px 
of view. From theoretical considerations, requirements regarding the molecular structure of 1g 
adhesive are deduced as regards its degree of polymerisation, branching, flexibility, and pres: 
of specific functional groups, which would ensure good wetting properties on the one ! 
and strong bonding on the other. The lack of literature references is regrettable in viev o. 
fundamental importance of this article. Wood technologists will find of interest the article by F 
E. C. Jahn, of the College of Forestry, Syracuse, N. Y., on the “Expanding Wood Utilisation thre 
the Use of Resins,’ in which the more important developments in wood-resin combinations 
discussed, such as glued wood construction, modified woods, paper overlaid lumber and plyw 
particle boards, paper-honeycomb sandwich material, various types of composite structures 
finishes. Particularly valuable as reference works are the articles by Drs. E. Ott andG. A. Heiber: 
of Food Machinery and Chemical Corp., New York, on ‘‘New Developments in Synthetic Re 
for Wood Adhesives,” where adhesives of the types of polyvinyl acetate, epoxy resins, phenol 


amine-aldehyde resins, and new resin adhesives and blends are reviewed. Detailed abstracts >f 1. 


recent literature references are appended and are a “must” for any worker in the field. 


The book is attractively printed on paper from a pulp blend containing Israeli Eucalyptus rostt 
produced at the American-Israeli Paper Mills, Hadera. Diagrams, photo-micrographs and electt 


photomicrographs complete the text where applicable and contribute to the general appearance o» 


the book. 
J. SCHMORA 
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